Scaling in Modulated Systems and Re-entrance of Order 
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We propose scaling invariance of d-dimensional modulated systems, in which modulation of order is 
produced by a frustrating long-ranged interaction decaying spatially with some power a of the inverse 
distance. Over the entire temperature range in which a modulation length is defined, measurable 
physical quantities acquire scaling exponents. One outcome of the scaling analysis is the existence, 
for a > d, of an anomalous state in which the tendency of matter to order increases with temperature. 
We suggest that this anomaly is responsible for a re-entrance phenomenon in two-dimensional dipolar 
frustrated magnetic systems. 
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the charge density in any type of strongly corre- 



Modulation of order is a general motive in chemistry, 
biology, and physics [J, 01 . The modulated order parame- 
ter^ay represent quantities as diverse as the spin density 

type of strc 

lated classical or quantum system 0, @, S Hi , the volume 
fraction of diblock copolymers, the concentration of am- 
phiphilic molecules and other chemical species 0, , or 
dipolar bosons in an optical lattice [lol | . However, modu- 
lated systems tend to show common characteristics such 
as the morphology of the various patterns and the occur- 
rence of transitions among them []|, as in liquid crystals 



[ill ] or two-dimensional melting phenomena V^. This 
tendency to common behavior indicates universal under- 
lying principles [tI. [sl. [Tst . 

In this Letter, we consider c?-dimensional systems 
where the spatial modulation of the order parameter is 
produced by the competition of a short-ranged interac- 
tion favoring local order and a weak but long-ranged frus- 
trating interaction decaying with some power a of the in- 
verse distance. We propose that these systems are scaling 
invariant along the entire temperature range where some 
characteristic modulation length L exists. Accordingly, 
measurable physical quantities are characterized by scal- 
ing exponents which only depend on fundamental sym- 
metries of the system, like d and a. One outcome of the 
scaling analysis is that modulation of order is only stable 
if a is less that an upper limit dependent on d. A further 
outcome is that, for a > d, the tendency of the mod- 
ulated system to order - measured by the compression 
modulus B - can increase with temperature. Among the 
systems that are predicted to show this anomaly are a 
model of two-dimensional Coulomb frustrated phase sep- 
aration Q and two-dimensional ferromagnets frustrated 
by the dipolar interaction 0, 0, @|- We report experi- 
mental data on perpendicularly ma gnet ized ultrathin Fe 
films on Cu(lOO) (a = 3 > d = 2) [liJil that confirm the 
anomalous growth of B with temperature. This growth 
of B goes along with a re-entrance [3l of the more or- 
dered stripe pattern 

Following Muratov 0, a universal effective Landau- 
Ginzburg- Wilson (LOW) Hamiltonian (per volume) Ti 
describing systems frustrated by a long-ranged interac- 



tion reads [1 
H\0[x 
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{x)(t>{x')Gx{\x-x'\)d^xd'^x' , (1) 



x being a vector in a D-dimensional space and Vd the 
volume of the system. Without the long-ranged interac- 
tion (A = 0), this functional has macroscopic phase sep- 
aration, with an order parameter 4>q appearing below a 
critical temperature Tc- The field 4'{x) is a scalar or- 
der parameter (rescaled with 0o)- The gradient term in 
the first line mimics the short-ranged exchange contribu- 
tion. It is inaccurate for sharp walls at low temperatures, 
where a more suitable expression exists [4] . For large sep- 
arations (|x — x'l), the kernel of the long-ranged interac- 
tion G\{\x—x'\) is assumed to behave asymptotically like 
jx — x'l^". a = l corresponds to the Coulomb interaction 
and gives rise to the Coulomb Frustrated Ising Ferro- 
magnet (CFIF) [6]. The kernel of the Dipolar Frustrated 
Ising Ferromagnet (DFIF) [23| has a more complicated 
structure but can also be treated with Eq. ([T]). We as- 
sume S> 1 and hence a mono-dimensional modulation 
of length L ^ 1 along a coordinate xi in a system with 
a large extent in d dimensions. In the remaining D ^d 
dimensions, the system is assumed to have a finite thick- 
ness 0{5) < L. 

At a given temperature T, we expect that a minor 
change of the modulation length L will have a limited 
effect on the order-parameter profile <j){x). We therefore 
postulate the invariance of (j){x) with respect to a scaling 
of all lengths by a constant u: 



{u-xi,u- L) = (j){xi,L) 



(2) 



i.e. (/) to be a function of one single variable Xi = we 
thus introduce (^) = (/)(a:i,L). Considering only the 
asymptotic behavior of Gx, we can now transfer the in- 
dependence of H to prefactors with suitable scaling 
exponents A. Per unit volume, the part due to the ex- 
change interaction scales like 0{J ■cji^-L^-') and the part 
due to the long-ranged interaction like 0(A-(/)g-L^^), with 
A,/ = — 1 (sharp walls between adjacent domains) or —2 
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(extended walls) and Aa = d — a. Some results of this 
scaling hypothesis are summarized in Table |T] and will be 
discussed in the following. 

- The generalized Hellmann-Feynman theorem [isj 
equates the derivative of the Helmholtz free energy per 
volume with respect to the parameter L with the canoni- 
cal average over all configurations (t>{xi) of the derivative 
of H with respect to L: 



dT{L, T) 
dl 



d_ 

dL 



(3) 
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The characteristic modulation length L\{T) is defined as 
the solution of the equation ^"^g^'"^^ =0 and, when finite, 
scales like (J/A)^^a with A^^ = 1/(Aa-Aj) (third and 
fourth columns). Note that, in general, walls are very 
thin at low temperatures (Aj = — f , third column) and 
broaden when the temperature is increased (A,/ = —2, 
fourth column) [i^] ■ This produces the cross-over of the 
exponent Aj;,^ observed in the third and fourth columns. 
Accordingly, scaling predicts a decrease of the equilib- 
rium modulation length with temperature. Particularly 
striking is the cross-over (and the corresponding decrease 
of La(T)) for the (£> = 3,rf = 2)-DFIF. For Aj = -1 
(T«0), we have A^^ =oo, which suggests an exponential 
dependence on J/A (confirmed by detailed computations 
0, S|; see Ref. [l9[ for d~ I and a = 2, which also give 
Alj^ =oo). At higher temperatures, L\ only depends lin- 
early on J/A. For J/A» f , this means a decrease of the 
modulation length by several orders of magnitude (20| . 
- The compression modulus Bx{T) 0, Hlj (also called 
Young modulus [2§|) measures the energy cost associated 
with deviations from the equilibrium modulation length 
L\(T) and is given by 



BxiT) 



-2 d^FiL,T) 



(4) 



iA(T) 



which, in virtue of the Hellmann-Feynman theorem, is 
equivalent to 



1 f d 



T \dL 



l4,{x) 



(5) 



Positivity of B\ is a necessary condition for the stability 
of the modulated state. The scaling analysis, applied 
J^n-iw. „. A(A/J)^I^ and 



to Eq. (O, shows that 



2-2A, 



T{{^'Hf)\Lx - (AVr)(A/J)^I^, so that positivity 



of Bx in the limit J/ A » 1 requires 

a - d < -Aj 



(6) 



This result generalizes the one obtained in Ref. [13] for 
T = and d = 1. When the modulated state is stable, 
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{D,d) 
DFIF(3,3) 
DFIF(3,2) 
CFIF(3,3) 
CFIF(3,2) 
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d — a 
1 

-1 
2 
1 



Bx ~ X-(j)Q- ^ (for A Bx i see fifth column). Note that 



TABLE I: Scaling exponents A^^ and A_g^ for general d and 
a and for some special cases. The asymptotic long-range in- 
teraction of the (D = 3, d = 3)-DFIF does not follow a simple 
\x — x'|~°-law. It therefore cannot be assigned a value of a 
and the exponents have to be determined separately. 



an anomalous L-dependence appears for d—a < (cf. the 
{D — 3,d — 2)-DFIF). Bx increases (by several orders of 
magnitude in particular cases) with decreasing Lx- This 
class of matter exhibits a state of potentially anomalous 
compression modulus, in which it becomes more difficult 
for the system to accommodate disorder involving devi- 
ations from the equilibrium modulation length when the 
temperature is increased; the system becomes less re- 
sponsive to fluctuations at higher temperatures. 

The validity of Eq. ^ is corroborated by a number of 
findings. 

- First, we have computed equilibrium profiles (j)cq{x, L) 
for the {d = 2)-DFIF, within the Mean-Field Approx- 
imation (MFA) of the functional Eq. ([T|) [13, Hi. At 
a given temperature, the MFA finds profiles that mini- 
mize the functional at a given L (black and light curves 
in Fig. 1(a)). Lx{T) corresponds to the profile (black 
curve) that minimizes 7i[(j)eqix , L) , T] with respect to L 
and is found to decrease with temperature (not shown 
[iot). All profiles virtually collapse onto one single "scal- 
ing profile" in Fig. 1(b) when plotted as a function of 

Figure 1(c) shows that the numerical Bq_{T) (black 
curve) for the (d = 2)-DFIF virtually coincides with the 
scaling prediction (light curve). In addition, it shows 
the anomalous growth with T predicted by scaling. Fig- 
ure 1(d) shows the almost perfect coincidence between 
the numerical and the scaling result for Bq{T) in the 
(d = 3)-DFIF and, in the inset, the (d = 2,3)-CFIF 
[a — 1). Bq{T) decreases smoothly with temperature 
and does not have any anomaly. 

- Second, we point out that Eq. ^ is trivially fulfilled 
at T = 0, as the order parameter field consists of seg- 
ments with length L and alternating values ± 1 . We have 
verified the stability of the scaling hypothesis and of the 
scaling dimensions at finite temperatures by introduc- 
ing small deviations from the strict mono-dimensional 
order, modifying (/)(xi) to (f){xi — ehkC0s{kx2) cos{hxi)). 
This replacement establishes an elementary excitation 
displacing domain walls along the a;i-coordinate while 
bending them along the a;2-direction and accounts for 
the experimental fact that stripes are never perfectly 
straight in real systems [l^] at finite temperatures. Such 
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FIG. 1: (a):(color online). Profile of the order parameter in 
a (d = 2)-DFIF {S = 1, fi/J = 0.08) for different values of 
L (0.7Ln, Lq and l.SLn). (b): Collapse of profiles onto one 
single curve in plot versus ^. (c): Numerical Bn{T) for the 
{d = 2)-DFIF (black curve, 5 = 1, Q./J = 0.08) together with 
the scaling result (light line), (d): Numerical and scaling 
curves Bn(T) for the (d = 3)-DFIF {S = 200, Q,/J = 0.08) 
and the [d = 2, 3)-CFIF (inset, 5 = 1,200, Q/J = 0.0001). 
Note that at low temperatures the wall is very thin and the 
numerical results are unstable. Light circles on the vertical 
axis indicate the values of B\ at T = obtained analytically. 

excitations are responsible e.g. for the Landau-Peierls 
instability that prevents order in D = 3 for d — 1,2 
0, S I3l- Averaging the perturbed functional H[(t){xi — 
thk cos(fca;2) •cos(/ia;i)), L, T] expanded up to order {ehkf 
and powers (h ■ k)" with n < 4 using the equipartition 
theorem (e^fe)^ w jq^h'^-^Xq-^-iki M a-^d summing over 
{h,k) < 1/L, we recover an effective mono-dimensional 
functional in which the exchange energy and the energy 
of the long-ranged interaction have, to leading order, the 
same scaling behavior, albeit with modified coupling con- 
stants J^J+0{T) and \^\+0{\'^). This shows that 
scaling is conserved at finite temperatures and that the 
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FIG. 2: Experimental Bq{T) for an ultrathin Fe film on 
Cu(OOl) obtained by dividing <?!)o(T) by Ln{T). The solid 
line results from fits to (j>o{T) and Ln{T) respectively and is a 
guide to the eye. The temperature range with stripe domains 
is shaded, the dark shade indicating mobile stripes. The two 
insets show SEMPA (Scanning Electron Microscopy with Po- 
larization Analysis) images of the labyrinthine and the stripe 
pattern. The length of the white bar is 10 /im, the thickness 
of the film 1.95 ML. 

scaling dimensions are not changed by small perturba- 
tions. In particular, the Landau-Peierls instability de- 
stroys long-range positional stripe order but the loss of 
order proceeds in such a way that the modulation length 
and the compression modulus Bx remain well-defined 
quantities also at finite temperatures [25J. 

- Third, the compression modulus Bq of perpendic- 
ularly magnetized ultrathin Fe films on Cu(lOO) 
shows the anomalous enhancement with temperature 
predicted by scaling. In these systems, a strong mag- 
netic anisotropy aligns the spins perpendicular to the 
film plane. Accordingly, the exchange interaction is frus- 
trated by the (isotropic part) of the long-range dipolar 
interaction, i.e. it is a (d = 2)-DFIF and has a character- 
istic modulation length [l], [3] , up to a temperature Tc 
at which 0o vanishes. The data points in Fig. 2 display 
Bq(T) as determined experimentally using the scaling 
result. The solid curve is obtained from fits to (f)o{T) 
and Ln{T). The decrease of 4>o{T) and the appearance 
of mobile stripes just below Tc makes the experimental 
determination of Bn{T) difficult. Note that in zero ap- 
plied magnetic field, the low-temperature state is striped, 
but transforms upon temperature increase into a more 
disordered labyrinthine pattern, displayed left in Fig. 2 
[l], This transformation has also been reported in 

a micrometer-thick film, a (rf = 3)-DFIF [l|. What is 
new in the ultrathin limit is that the labyrinthine pat- 
tern transforms back into the more ordered stripe pat- 
tern (pictured right in Fig. 2) when the temperature is 
increased. This re-entrance of order proceeds in the tem- 
perature region (light shaded area) where the widespread 
local configurations in the labyrinthine pattern with devi- 
ations from the equilibrium stripe width are increasingly 
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penalized by the increasing Bq{T). This suggests that 
the unusual re-entrance of stripe order is a response of 
the system to the extra elastic energy introduced by the 
increase of Bq (T) . Note that the ex per imental tempera- 
ture dependence of the stripe width [15| is smooth in the 
temperature range where the labyrinth-to-stripe transfor- 
mation proceeds, which confirms the predicted insensitiv- 
ity of Lfi to deviations from mono-dimensional order. In 
the vicinity of Tc (dark shaded area), the stripes become 
mobile [15[. This mobility coincides with the range of 
temperatures where Bq{T) is expected to decrease and, 
accordingly, can no longer counteract thermal excitations 
which destroy stripe order. 

In conclusion, we have presented a scaling procedure 
that is generally applicable to any pattern-forming sys- 
tem and allows re- writing the Hamiltonian ([1]) as a func- 
tional of a scaling-invariant scalar field and rescaled cou- 
pling constants J i-^ JL^' andAi-^A-L^^. The tempera- 
ture dependence of La, which may vary by several orders 
of magnitude, is thus expected to propagate in a non- 
trivial way to any thermodynamic observable. Among 
other possible implications, scaling predicts a yet unde- 
tected state of matter with a compression modulus in- 
creasing with temperature. This anomaly provides a sim- 
ple explanation for an anomalous sequence of patterns 
observed in ultrathin ferromagnetic films. We show that 
a necessary condition for this anomaly is a > d, which is 
fulfilled by the ultrathin limit of dipolar frustrated Ising 
fcrromagnets (or equivalent models, see Ref. but 
not by dipolar frustrated films that arc not strictly two- 
dimensional or in Coulomb frustrated systems. This ex- 
plains why, despite the huge amount of experimental and 
theoretical work, no other instances of re-entrance have 
ever been observed in the many other pattern-forming 
systems [l|. So far, predictions and studies of topological 
transitions although extensive, have neglected possi- 
ble anomalous behavior of the compression modulus. Our 
analysis indicates the need to refine the various scenarios 
for topological phase transitions. In addition, we antici- 
pate that tuning B\ with an appropriate choice of tem- 
perature, thickness - and probably applied magnetic field 
- might provide new ways of producing two-dimensional 
melting transitions and allow investigating their elusive 
topological character as well as the controlled cross-over 
to more conventional three-dimensional melting. We also 
anticipate that other re-entrance phenomena - e.g. in 
an applied magnetic field - might arise because of the 
anomalous compression modulus. 
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